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Abstract 

According to previous work, topological defects expand exponentially 
without an end if the vacuum expectation value of the Higgs field is of the 
order of the Planck mass. We extend the study of inflating topological defects 
to the Brans-Dicke gravity. With the help of numerical simulation we inves- 
tigate the dynamics and spacetime structure of a global monopole. Contrary 
to the case of the Einstein gravity any inflating monopole eventually shrinks 
and takes a stable configuration. We also discuss cosmological constraints on 
the model parameters. 
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I. INTRODUCTION 



For the last decade spacetime solutions of gravitating monopoles has been intensively 
studied in the literature |l|-§. This originated from rather mathematical interest in static 
monopole solutions. In both cases of global monopoles [jl],|| and of magnetic monopoles || , 
static regular solutions are nonexistent if the vacuum expectation value (VEV) of the Higgs 
field 7] is larger than a critical value r] sta , which is of the order of the Planck mass, mp\. The 
properties of monopoles for 77 > ?7 sta had been a puzzle for some time. 

In connection with the above issue, it has been claimed by several authors independently 
that monopoles expand exponentially if 77 > 0(rrip\) 0-0. Among them, Guendelman 
and Rabinowitz j3J assumed the simplified model of a global monopole under the thin-wall 
approximation, while Linde || and Vilenkin |J made qualitative arguments on any kinds 
of topological defects and argued that their core would inflate if the VEV of the Higgs 
field is large. Later the full evolution equations for monopoles were numerically solved 
I?]!!, supporting the previous arguments as a whole. Furthermore, those numerical analyses 
revealed an unexpected result that the critical value, rj in f, above which inflation occurs in 
the core, is larger than ?7 sta , and stable but nonstatic solutions exist for ?7 sta < r\ < r] in{ [Q. 
For example, in the case of global monopoles, we find i] sta = 0.20m P1 and i] in f = 0.33mpi ||. 
Global spacetime structure of an inflating monopole has also been discussed in . 

In this paper, we extend the study of monopole inflation to the Brans-Dicke (BD) theory. 
If the BD field exists, the inflationary universe may not expand exponentially but expand 
with a power law even in the presence of an effective cosmological constant |10|JTl[|. For 



example, in extended inflation [jnj, this slower expansion was expected to solve the graceful 
exit problem of old inflation [12]. We thus expect that the BD field also affects the dynamics 
and global spacetime structure of inflating monopoles. 

We would also like to discuss whether the present model can be a realistic cosmological 
model or not. Once inflation happens, the spatial gradients of the fields become negligibly 
small on the scale of the observed universe. Therefore, just as other inflationary models, we 
have only to care about quantum fluctuations. This issue resolves itself into the commonplace 
analysis of density perturbations in the inflationary universe. Here we adopt the result 
obtained by Starobinsky and Yokoyama [[D|], who derived the general expression of the 
density fluctuations in the BD theory, taking account of the isocurvature mode as well as 
the adiabatic mode. 

The rest of the paper is composed of two independent analyses. In the first part (§2), we 
investigate the dynamics and global spacetime structure of an inflating global monopole. In 
the second part (§3), we discuss cosmological constraints on this model, which come mainly 
from density perturbations. We use the units c — h—1 throughout the paper. 

II. DYNAMICS AND SPACETIME STRUCTURE OF A GLOBAL MONOPOLE 

The Brans-Dicke-Higgs system, which we consider here, is described by the action 

<E> uj ,„ ^9 1 



(2.1) 



with 



2 



V(V) = ^ 2 -r/ 2 ) 2 , *e\/W^, (2.2) 

where $ and \l/ a are the BD field and the real triplet Higgs field, respectively, ui and A are 
the BD parameter and the Higgs self-coupling constant, respectively. 

Let us begin with a discussion of the fate of an inflating topological defect. Once infla- 
tion begins, the core region can be approximated by the flat Friedmann-Robertson- Walker 
spacetime: 

ds 2 = -dt 2 + a(t) 2 dx 2 , (2.3) 



which yields the solution of extended inflation fT^JTOjl 



a(t) oc ( 1 + —J , $(t) oc f 1 + -M , (2.4) 



where Hi = J 8TcV(0)/3§i is the Hubble parameter at the beginning of inflation, and a 2 = 
(2uj + 3)(6u + 5)/12. Hence, the effective Planck mass, 

i 

m P i($) = oc a^+i , (2.5) 

continues to increase until inflation ends. This implies that, even if ?y/mpi($) is large enough 
to start inflation initially, it eventually becomes smaller than the critical value (= 0.33) 
We thus speculate that any defect eventually shrinks after inflation. 

In order to ascertain the above argument, we carry out numerical analysis for spherical 
global monopoles. The coordinate system we adopt is 

ds 2 = -dt 2 + A 2 (t, r)dr 2 + B 2 (t, r)r 2 (d6 2 + sin 2 6d V 2 ). (2.6) 

For the matter fields, we adopt the hedgehog ansatz: 

M> a = M>(t ;r )f a , f a = (sin0cos<^,sin0sin<£,cos6>). (2.7) 

For the initial configurations of the Higgs fields and the BD field, we assume 

= 0,r) =?7tanh(^) with 5 = ^2-, *(t=0,r) = 0, (2.8) 

$(t = 0, r) = $j = const., $(t = 0, r) = 0, (2.9) 

where an overdot denotes d/dt. 

Although the numerical code developed by one of us || worked well for monopoles in the 
Einstein theory, for the present system which includes the BD field it sometimes does not 
keep good accuracy for sufficient time. We therefore improve our numerical method, mainly 



based on the idea of Nakamura et al. [15| The basic equations and our improved numerical 
method are presented in Appendix. 

To illustrate our results, we take u — 1, A = 0.1 and r//mpi($j) = 1. The evolutions of 
^ and 7//mpi($) are shown in Figure 1. Clearly the above arguments are verified by this 
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numerical result: as r}/mp\($) gets close to the critical value, the monopole stops expanding 
and turns to shrink. In Fig. 2 we plot the trajectories of the position of \P = T]/2 for several 
values of to. The curves indicate that a monopole do not shrink toward the origin but tends 
to take a stable configuration. 

Next, we shall examine the metric outside the monopole. Before we consider monopoles 
in the BD theory, let us review previous work on the outer solution in the Einstein theory 
briefly. The approximate expression for the metric outside a static monopole was found by 
Barriola and Vilenkin (BV) 

ds 2 = - ( 1 - A - dT 2 + ( 1 - A - 1 dR 2 + R 2 (d6 2 + sin 2 Bdy 2 ). (2.10) 

\ m-p\R J \ m p \R J 

where M c is an integration constant and A = 87T7/ 2 / 'm P \. Later Harari and Lousto pointed 
out that the mass parameter M c is negative and therefore a global monopole has a repulsive 
nature ||. For large R, the metric ( |2.10| ) can be approximated by 

ds 2 = - (1 - A) dT 2 + (1 - A)" 1 dR 2 + R 2 {d6 2 + sin 2 9d(p 2 ), (2.11) 



and therefore Air A is interpreted as a deficit solid angle. Although the expression ( p. 11 ) looks 



singular when A = 1, Cho & Vilenkin || showed that it is just an apparent singularity and 
( [2. 11| ) also applies to the case A > 1. 

Now we would like to discuss the properties of the outer spacetime by comparing it with 
the BV solution ( |2.10|) . As a coordinate-independent variable, we adopt the Misner-Sharp 



mass | 16| , which is defined as 



M^^il-g^R^ (2.12) 

where R = ^fgee is the circumferential radius of the spacetime. Note that M has a dimension 
of length, and m P \M is a mass in a usual sense. The BV solution is characterized by 
M = M C + AR with M c = M c /m P \. 

Figure 3 shows a plot of M/R. In the case of the Einstein theory (a), M/R seems 
to converge into A/2 (=Att in this case) at large R, confirming that the metric is well 
approximated by ( |2.10| ). The result with the BD gravity looks similar, but a different 
character is that M/R decreases with time Gvcn sit ct fcir region. This behavior is consistent 
with the BV solution in the Einstein theory: if we take account of the time-variation of 
m-pi in the BV expression ( |2.10| ), a deficit solid angle effectively decreases and ( P-IOQ also 



approximates the solution in the BD theory. 

Our results obtained so far are rephrased as follows: our numerical integration of the evo- 
lution equations verifies as a whole the arguments with considering only the time-variation 
of mpi($). That is, the metric in the core of a monopole is well approximated by ( |2.4D , 
and the metric outside the core by ( |2.10|) with changing mpi. Although we have chosen 



small to to illustrate the effect of the BD field clearly, if we take u > 500 as constrained by 



observation ||17|| , the spatial gradient of $ is much smaller, justifying more the arguments 
with neglecting the spatial configuration of $. 

The result that inflation eventually ends stems merely from the change of the local values 
of mpi($). We may therefore extend this result to models of other topological defects. 
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III. COSMOLOGICAL CONSTRAINTS 



In this section we discuss constraints for this model to be a realistic cosmological model. 
As we mentioned in §1, the main constraint is from density perturbations. In the present 
model matter fields are composed of the BD field $ and the Higgs field \I/ a . For the Higgs 
field, we only calculate fluctuations in the radial direction 5^, because fluctuations along 
the other directions do not contribute to growing modes. 

Here we analyze the field equations in the Einstein frame by use of a conformal trans- 
formation, following previous work |ll|,13|. A wide class of generalized Einstein theories is 



described by the Einstein-Hilbert action with scalar fields: 



S 



d x 



2^ - 2 W) 



e —— (V#) 2 



-P**V(V) - U((j)) 



with k 



87rG, 
(3.1) 



where f3 and 7 are constants, and <fi is the redefined scalar field. 

In the BD theory, on which we concentrate in this section, the original action in the 
Jordan frame (12. II) is transformed into (13.11) with 
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u{<p) = 



2 V 2cu + 3 ' 

through the conformal transformation and the redefinition of the scalar field 

9p,v = -j-9nv> VlK<j) = In — , 



(3.2) 



(3.3) 



with $0 being a constant. From 



the Planck mass in the Jordan frame is written as 



m P i($) 



$ e 



7«</>/2 



(3.4) 



We set = at the present time, t = t , identifying a constant y^o with the Planck mass at 
present. Because the evolution of the BD field after inflation is negligibly small for 7 2 1, 
we may identify $ at the end of inflation, <3>y, with $ (ie., <frf = <f)Q = 0). In the following 
we sometimes use the symbol mpi($) instead of 0, and define mpi j0 = mpi($ ) = mpi($/). 

Hereafter, we omit a hat, which has denoted quantities in the Einstein frame. Taking 
the background as the spatially flat Friedmann- Robertson- Walker spacetime (|2.3|) , the field 
equations for the homogeneous parts read 



H 2 



a\ k' 
a) ~ Y \2 6 



1 



(3.5) 



dV 



^ + m^> - + e~ 1K<p — = 0, 

d^ 



+ 3H<j) + ^V 7K< ^ 2 - 2 7 Ke- 2W 1/ = 0. 



(3.6) 
(3.7) 
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The conditions for slow-roll inflation <C 1 3 J^r\I/ 1 , \<p\ <C \3H<p\, etc.) are equivalent to 



Y « J, 



m P i($)y^ 



V 



< V/487T, 



m P i($) 2 V^ 



< 24vT, 



(3.8) 



where = ci^/d\l/. For the potential (|2.2| 
in the range < ^ < rj are equivalent to 



the second and the third inequalities in ( |3.8| 



e = 



67rr/ 



m P i($) 



> f, ^ < r? 



e 

f - - 

2 



(3.9) 



where terms of higher-order of e have been neglected in the second inequality. 

One of the distinguished features of this model is that there are two scenarios of exiting 
from an inflationary phase and entering a reheating phase. This is illustrated by the slow-roll 
conditions ( |3.9|) . When the condition e _1 ^> 1 (a more precise condition is rj/mpi(&) > 0.33) 
breaks down, inflation stops globally. Before this time, many local regions with the present- 
horizon size enters a reheating phase when the second condition \1/ -C breaks down. The 
second scenario is just like that of standard topological inflation or other slow-roll inflationary 
models. In the first scenario a microscopic monopole might remain in the observable universe. 
Unfortunately, however, this possibility turns out to be ruled out because the spectrum of 
density fluctuation is tilted excessively if r\ is too close to the critical value 0.33m P1 |TJ) and 
because r)/mpi($f) must be larger than 1 from the COBE normalization as will be seen later 
in Fig. 4. Thus we only consider the second (standard) reheating scenario below. 

When inequalities (|3.9|) are satisfied, the field equations are approximated by 



H 2 



2-jK(f>y 



(3.10) 



3H^ = -e 



~"fK<f) 



dV 



(3.11) 



which yield solutions , ], 



3H<P = 2 7 Ke- 2W V, 



V a > 2 7 ' 



(3.12) 



(3.13) 



N = = k 2 

2 7 2 hf V> 



7T 



m 



P1,0 



2r? 2 In -4- + ^ 2 - 



(3.14) 



Starobinsky and Yokoyama [JT3 1 derived the amplitude of density perturbation on comov- 
ing scale I = 2ti fk in terms of Bardeen's variable $a [19] as 



48V 



25m P1 ($) z 



g7r y2 , e 2^N _ 1)2 

+ 



4 7 2 



(3.15) 
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where all quantities are defined at the time tk when /c-mode leaves the Hubble horizon, i.e., 
when k = aH . Since the large-angular-scale anisotropy of the microwave background due 
to the Sachs- Wolfe effect is given by 5T/T = $,4/3, we can constrain the values of A and 
77/rapi by the COBE-DMR data normalization f2C§, 



i$ A = 1(T 5 , (3.16) 

on the relevant scale. To calculate ( p,16|) explicitly, we have to relate scales during inflation 
and present scales. The number of e-folds iV between t = tk and t — tf is given by 



iV fc = (l-7 



2\-l 



68-hi— - + ln^ + 21n^V--ln^- | , (3.17) 
^ a H m Pli0 V 4 3 p ^ 



where p r h is the energy density when reheating completes, which is not determined without 
specifying a reheating model. Here we choose Nk =ao H = 65 typically. Then the correspond- 
ing values of <ft an d ^ are determined from ( |3.13j ) and ( |3.14| ). Figure 4 shows the allowed 



values of A and 77/mpi for 7 = 0.045 (00 = 500). The concordant values are represented by 
two curves, which is a distinguished feature for the double- well potential (|2.2|). Unfortu- 
nately, fine-tuning of A ( < 10~ 13 ) is needed, just like other models fpL3|| . The condition of 
A < 10 -13 in the present potential was also found in [pi]. The constraints for uo > 500 are 



practically no different from those in the Einstein gravity. 

Using the relation dink = da/a + dH/H and ( |3.10| )-( |3.12| ), we obtain 

n 1 = rfln$j ^ 3m P1 ($)% 2 | m P1 ($) 2 ^ 3 (3 lg) 

d\nk 8ttV 2 AtcV 

The spectral indices are plotted in Fig. 5; the two lines correspond to the two lines of 
COBE-normalized amplitudes in Fig. 4. As ( |3.18| ) suggests, the deviation from the Einstein 
theory is only 67 2 = 0.01 for u = 500. Therefore, relatively large shift from n — 1 is caused 
not by the BD field but by the double well potential (|2.2|). 

Besides the amplitude of density perturbation, we have to check other conditions for 
successful inflation. First, if we take quantum fluctuations S^/q into account, the classical 
dynamics is meaningless unless 

U > Wq = __. (3.19) 

We check the value of ^ at k = aoHo which satisfies the COBE normalization, finding that 
( p,19| ) is always satisfied. 

The second is the condition that classical description of spacetime is valid. We examine 
this because we have considered the case of r\ = mpi, which may be critical. We require 

V(%) < m P1 ($,) 4 . (3.20) 

This condition is also satisfied for rj = mp^o and A < 10~ 13 . In the end the two conditions 
do not add new constraints on the model. 
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Now, let us discuss the detectability of relic monopoles in this model. First, we investi- 
gate how the long-range term of the monopole mass density, p gm (R) = r] 2 /R 2 , contributes 
to the cosmological mass density, where R is the distance from a monopole core to an ob- 




server. Hiscock |22| derived the condition that the averaged mass density of multiple global 
monopoles does not exceed the cosmological mass density, i.e., (p gm ) < Po- His result is 
equivalent to 

AT i /™ \ 

(3.21) 

where N/Vh is the number of monopoles in the the present horizon volume, Vh = 
(47r/3)-£f(7 3 , and VLq is the present density parameter. This indicates that the existence 
of even one monopole is critical for rj = mpifi- To see it more closely, we consider the condi- 
tion that the local density by one monopole does not exceed the cosmological mass density, 
i.e., p gm (R) < po- This leads 

R>**jjl(-JL-), (3.22) 



'Q \mpi^ / 

which shows that the center of a monopole should be located farther than the horizon. A 
more stringent constraint is obtained by considering the difference of p gm at both ends of 
the observed universe: 



rf_ r, 2 _ Ar] 2 H 1 

-Ho 

The condition Ap gm / p < 1CT 5 leads 



APsm " (R - Ho 1 ) 2 (R + Ho 1 ) 2 ~ R* ' 1 J 



R> 1 -^(^-Y. (3.24) 

We now argue the plausible value of R realized in this model. By assumption our universe 
was once in the core of a monopole with \1/ = and the Higgs field started classical evolution 
only after it had acquired an amplitude $ > c^q. The number of e-folds after this epoch 
is approximately given by 



N = 2*[-H-) In » . (3.25) 




For 7] = m-p\fi and A = 10~ 13 , we find iV = 103. That is, the initial core radius becomes 
e i03-65 ~ 10 16 times as large as the present horizon. We therefore conclude that the core 
of a monopole is located so far from the observed region that even its long-range term 
Pgm — rj 2 / R 2 does not exert effective astrophysical influence. Because this conclusion is 
irrelevant to the long-range term, it suggests that a local monopole in the same model is not 
detectable either. 
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IV. SUMMARY AND DISCUSSION 



We have considered inflating monopoles in the BD theory. First, we have investigated the 
dynamics and spacetime structure of a global monopole. In the Einstein theory, monopoles 
undergo eternal inflation if and only if r)/mp\ > 0.33. In the BD theory, on the other hand, 
the ratio i]/mp\(§) decreases until inflation ends. That is, if r)/mp\{$) is large enough 
initially, a monopole starts to expand; however, after the ratio becomes smaller than the 
critical value (= 0.33), the monopole turns to shrink. This behavior is quite different from 
that in the Einstein gravity and gives another mechanism to terminate inflation, which, 
however, is irrelevant to our universe. 

Secondly, we have examined cosmological constraints on the parameters, which stems 
mainly from density perturbations. The COBE normalization requires r]/mp\fi > 1, 
A < 10~ 13 for uo > 500, implying that inflation in our universe must have ended in the 
standard manner as \l/ reaches r]. We also calculate the spectral index, which turns out to 
be typically around n = 0.9. 

Although we have concentrated on the BD theory, finally we make a brief discussion on 
the behavior in general theories ( |3.1| ). We still assume U(<f)) = 0. The slow-roll conditions 
for the models are 



£ and \^\ « 3, 



K 



V 



<C Vo and 



\ 



6/3 



7 



(7^0), 



V 



V 



<3 



(4.1) 



and the behavior of cf) is given by 



(4.2) 



These show that, if (3^ > 0, 70 increases and the slow-roll conditions eventually break down. 
We may therefore conclude that the finite duration of topological inflation is a general nature 
in generalized Einstein theories. 
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APPENDIX A: FIELD EQUATIONS AND NUMERICAL METHOD 



In this Appendix we explain how we solve the field equations for a global monopole. As 
we shall show below, we improve our previous schemes pj so as to keep better accuracy. 



The variation of fl2.1| ) with respect to g^, $ a , and $ yield the field equations: 



G 



ft]' 



K 



T 

dV(V) 
2uj + 3'' 



with 



1 



(Al) 
(A2) 
(A3) 

(A4) 



In the following we shall write down the field equations (|A"l|), (^2), and ( |A3|) with the 
metric ( |2.6|) and the hedgehog ansatz (2/7 ). To begin with, we introduce the extrinsic 



curvature tensor K. 



Kl 



4. *S 



A K 
B 1 



Kl 



(A5) 



Next, following Nakamura et al. flgl , we define the following variables to guarantee the 
regularity conditions at the center: 



A-B 



k 



Kl - Kl 



(A6) 



and introduce a new space variable, x 



Because only \l/ is an odd function of r due to 



the hedgehog ansatz, we adopt a variable ip = ^f/r instead of ty. 

There are several advantages in the above procedure of Nakamura et al. |15| The first 
is that any diverging factor at r = like 1/r does not appear in the basic equations below; 
hence we do not have to treat the r = point separately. The second is that the new 
variables a and k are replaced with differences between comparable quantities like A — B at 
r « 0, which would generate numerical errors due to a finite decimal of computers. The third 
is that a derivative at r w is approximated by a finite difference much more accurately. 
To see this, let us imagine a function F(r) which is expanded as 



F(r) = c + c\r + c 2 r 2 + c 3 r 3 + 



(A7) 



around r = 0. If we take the central difference of F with respect to r, its derivative at 
r = Ar is approximated as 



Ar1 



F(2Ar) - F(0) 
2Ar 



0(Ar 2 



(A8) 



that is, it is approximated only up to the c 2 -term. If F is given as an even or odd function 
of r and if we use a space variable x, on the other hand, we can expand it as 
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± even 

F dd 



c + c 2 x + C4X + ■ ■ 
C\ + c 3 x + c 5 x 2 + • • 



(A9) 
(A10) 



If we take the central difference with respect to x at x = (Ar) 2 , for an even (odd) function, 
it approximates dF/dx up to c 4 (c 5 ), and moreover guarantees c 5 = (c 6 = 0). 

Furthermore, we define the following auxiliary variables to make the equations have a 
first-order form: 



B' 



n 



(All) 

(A12) 
(A13) 



where a prime and an overdot d/dx and <9/<9t, respectively. A full set of dynamical variables 
is a, B, C , ip, £, $, 5, K, k, zu, and II. The time derivatives the first seven variables are 
given by the definitions above: 



a = -aK r r + Bk, 
B = —BKq, C = K b e 
ip = zu, £ = zu', 

$ = n, s = n'. 



(A14) 
(A15) 
(A16) 
(A17) 



Note that the value of Kg' in ( |A15| ) can be determined by the momentum constraint (|A19|) , 
which gives a more accurate value than a finite difference of Kg. 
Now we write down the field equations (|A~T|) , ( |A2| ) and (|A3| ) as 



-G\ = 



K 2 - k 2 x 



2^,2 



ix 



A'C 3C 2 



A 



4 



A' a(A + B) 



87T 



xzu 2 2x£ . _ ,. ijj 2 ijj 2 



2 

n full 



Gtr 

4r 



Kf + k 



A 2 
k) + 



4 



xC 



A 2 $ 
2ttzu 



2A 2 B 2 
~ ( ^ 



+ xi. 



2$ 



A' 
~A 



2C 



$ 2$ 



n< + S (^ + AT) 



n\ = K 



K 2 + 2x 2 k 2 8 



7T 



V) 



12ttT* 
(2cj + 3)$ 
4 



A 2 $ 



'A' 



(A18) 
(A19) 

(A20) 



n r r -n b 6 



8vr 



A 2 $ 



k — kK + 



A 3 



2AC" - 2A'C + aC + a 



a 

B 2 



2 ,4 
~2 



B 



ip 2 a(A + B) 
B 2 



kU 4 
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w = Kuj + 



n = kti 



A 2 
4 

A 2 



A 2 



with 



xzu 



xz. 



A 2 



2u; + 3' 



4V. 



(A22) 
(A23) 

(A24) 



In order to set up initial data, we assume A(£ = 0,r) = B(t = 0,r) = 1 besides the 
matter configurations ( |2.8| ) and (|2.9|), and solve the constraint equations ( |A18| ) and ( |A19| ) to 
determine K(t = 0, r) and k(t = 0, r). Equation (|A20|) , ( |A2T| ), Q52g ), and (|A2l ) provide the 
next time step of K, k, w, and II, respectively. The constraint equations ( |A18[ ) and ( |A19| ) 
remain unsolved during the evolution and are used for checking the numerical accuracy. 
Through all the calculations the errors are always less than a few percent. 
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FIG.l. Dynamics of a global monopole with u — 1, A = 0.1, and rj/mpi(^i) = 1. We plot 
^ and 77/mpi($) in (a) and in (b), respectively. Because we omit the far region where W is 
almost constant in (a), the scales of abscissas in (a) and in (b) are different. 
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FIG. 2. Dependence of the evolution of a global monopole on u. We set A = 0.1 and 
r)/mp\($i) = 1. We plot trajectories of the position of \& = r]/2. 
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FIG. 3. Behavior of the Misner-Sharp mass outside the core of a monopole. We take, for 
reference, the Einstein gravity in (a), and u — 1 in (b). 
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FIG. 4. Concordant values of A and 77/m.py) with CO BE- normalized amplitudes of density 
perturbations. 




FIG. 5. The spectral indices of density perturbations. These two index curves correspond 
to the two amplitude curves in Fig. 4. 
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